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Abstract 

^5 ■ We analyze time-of-flight absorption images obtained with dilute Bose-Einstein con- 

Q . densates released from shaken optical lattices, both theoretically and experimentally. 

^ O^ l We argue that weakly interacting, ultracold quantum gases in kilohertz-driven opti- 

cal potentials constitute equilibrium systems characterized by a steady-state distri- 
^N . bution of Floquet-state occupation numbers. Our experimental results consistently 

Q«^ ! indicate that a driven ultracold Bose gas tends to occupy a single Floquet state, just 

IjT) ] as it occupies a single energy eigenstate when there is no forcing. When the driving 

amplitude is sufficiently high, the Floquet state possessing the lowest mean energy 
does not necessarily coincide with the Floquet state connected to the ground state 
of the undriven system. We observe strongly driven Bose gases to condense into the 
^S| ' former state under such conditions, thus providing nontrivial examples of dressed 

matter waves. 
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1 Introduction 



There is a growing interest in ultracold atoms confined in time-periodically 
driven optical lattices. It was pointed out already in 1997 that a metal- 
insulator-like transition undergone by ultracold atoms in quasiperiodic optical 
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lattices c an be induced by changinK the amplitude of a time-periodic driv- 
ing force ( jDrese and Holthaua . 119971 ) . But the field became truly active only 
10 years later, after the unambiguous experimental observation of dynamical 
tunneling suppression with Bose-Einstein condensates, and even of driving- 
i nduced reversal of t he sign of the tunneling matrix element, by the Pisa group 
( JLignier et al.l . 120071 ). and after related single-par ticle experinients p erformed 



with driven double- well potentials in Heidelberg fJKierig et al.l . l2008l ). Several 



groups then have monitored the dynamics under the combined action of both a 
homogeneous time-indepen dent and a time-per i odic driving force, ei ther with 



non-condensed cold atoms (llvanov et al.l. l2008l: lAlberti et al.l. 20091 ). or with 



Bose-Einstein condensates (ISias et al.l . l2008l : iHaller et al.l . l2010l): such dynam 



ics ma y be viewed as an analog of photon-assisted tunneling (JEckardt et al. 



2005al ). The Pisa group subsequently also reported the coherent control of 
the superfiu id-to-Mott insulator transition in shaken three-dimensional opti- 
cal lattices ( jZenesini et al.l . |2009|) . 



The recent detection of spontaneous breaking of time- rever s al sym metry with 
fast-oscillating triangular optical lattices by IStruck et al.l (120 111 ) constitutes 
a further milestone in this line of research. Moreover, photon-assisted tun- 
neling in a strongly correlated Bose gas has very recently been investigated, 
poss ibly allowing for applications to topological physics and quantum comput- 
ing ( jMa et al.l . 1201 ll ). Also active control of correlated tunneling in ac-driven 
optical superlattices has been o bserved, providing a novel approach to the 
realization of XX Z spin models (IChen et al.l . 1201 ll ). 



On t he theoretical side, several st udies indicate promising future perspec- 
tives: lEckardt and HolthausI ( l2008bl ) have suggested to perform avoided-level- 
crossing spe ctroscopy with ultrac old Bose gases in strongly driven optical lat- 
tices, while lEckardt et al.l (|2010[ ) have worked out a proposal for simulating 
frustr ated quantum antiferro n aagne tism with oscillating triangular lattices. 
Later iTokuno and Giamarchil (120111 ) have pointed out that a small periodic 
phase modulation of an o ptical lattice can g ive direct access to the system's 
conductivity. Furthermore, iTsuji et al.l ( 1201 ll ) have argued that ac forcing may 
even change the interparticle interaction from repulsive to attractive, thereby 
enabeling one to simulate an effectively attractive Hubbard model with a tem- 
perature below the superconducting transition temperature. 



On a first superficial glance, ultracold atoms in shaken optical lattices, with 
shaking frequencies on the order of a few to tens of kilohertz, might appear 
as typical examples of nonequilibrium systems. But this naive view is not 
correct. Quantum systems governed by a time-periodic Hamiltonian possess 
a particular basis, consisting of Floquet states, with respect to which one 
finds occupation numbers which remain constant in time, with the periodic 
time-dependence already being incorporated into the basis states themselves. 
Thus, a weakly interacting quantum gas subjected to time-periodic forcing. 



and possibly exposed to some sort of noise, will be statistically characterized by 
a time-independent distribution of Floquet-state occupation numbers. What 
is that equilibrium distribution? 

In the present tutorial article we develop this line of thought, and discuss ex- 
perimental data, obtained with Bose-Einstein condensates in kilohertz-shaken 
optical lattices, which we interpret as evidence for Floquet-based dynamics. 
We proceed as follows: In Sec. [2] we review basic elements of the Floquet pic- 
ture, and speculate about the existence of analogs of the Bose-Einstein and the 
Fermi-Dirac distribution function for isolated, time-periodically driven quan- 
tum gases. We then turn in Sec. [3] to t he realization o f perio dically shaken 



optical lattices, extending the analysis of iMadison et al.l (jl997l ). and pay par- 
ticular attention to the distinction between the co-moving frame of reference, 
in which the trapped quantum gas experiences a spatially homogeneous iner- 
tial force, and the laboratory frame, in which the measurements are performed . 



In Sec , m we briefly introduce the driven Bose- Hubbard model (JEckardt et al. 



2005bl ). In Sec. [5] we calculate momentum distributions expected in time-of- 
fiight absorption imaging for matter waves occupying a single Floquet state. 
We consider two particular situations: An interacting gas periodically driven at 
high frequencies, and a noninteracting gas subjected to a force with arbitrary 
time-dependence. Here we also point out certain characteristic differences be- 
tween time-of-fiight absorption images obtained from shaken lattices utilizing 
the inertial force, and the corresponding images acquired when the lattice 
remains at rest, while the forc e is exerted, e.g., t hrough the harmonic modu- 



lation of a levitation gradient ( iHaller et al.l . 120101 ). In Sec. [6] we tie the various 



strands together and report our measurements, which suggest the existence of 
Floquet condensates occupying the "lowest" available Floquet state. While we 
consider only one comparatively simple experimental setting, it appears likely 
that our findings are of a more general nature. In the final Sec. [7] we briefly 
spell out our main conclusions. 



2 The quest for Floquet condensates 



The notion of Bose-Einste in condensation of an almost ideal gas (JEinstein 



1924J . Il925l : iLondoru . Il938[ ) appears to be inextricably linked to the energy 
eigenstates of some single-particle Hamiltonian: When the temperature of a 
gas of weakly interacting Bose particles becomes lower than a certain crit- 
ical temperature, the gas "condenses" into the single-particle ground state, 
meaning that this ground state becomes macroscopically occupied. It might 
seem that this concept cannot be applied to explicitly time- dependent sys- 
tems, for which, in general, there are no stationary states. Yet, there is one 
notable exception: If a single-particle Hamiltonian H^^\t) is periodic in time, 
so that H^^\t) = H^^^t + T) for some period T, then Floquet's theorem 



( jFloquetl . Il883l : iKuchmentl . Il993[ ) suggests the existence of a set of distin- 
guished solutions lipnit)), known as Floquet states, to the time-dependent 
Schrodinger equation ihdtliljit)) = H^^\t)\^{t)); these Floquet states are in 
many ways analogous t o the usual energy eigenstates of tirne-independent 



Hani i ltonian operators (jShirlev 



19731 : iBreuer and Holthausl . 11991 
ticular form 



1965; Zel'dovich. 1967; Ritus. 1967; Sambe 



Chu and Telnovl . |200J) . They have the par- 



iV^n(^)) = \un{t))exp{-ient/h) 



:i) 



where the functions \un(t)) = \un(t + T)) inherit the period T of the underly- 
ing Hamiltonian; the index n denotes a set of quantum numbers specifying the 
state. The quantities En, which determine the phase factors accompanying the 
time-evolution of the Floquet states ([1]) in the same manner as energy eigenval- 
ues of time-independent systems determine the phase factors ex pressing the 
evolution of their energy eigenstates, are dubbed quasienergies (IZeFdovichl . 
19671 : lRitusl . ll967h . 



There is a compelling reason for considering these Floquet states. Namely, 
at each instant to the set of Floquet functions {|Mn(^o))} is complete in the 
physical Hilbert space on which H^^\t) acts. Hence, any solution lipit)) to the 
time-dependent Schrodinger equation admits an expansion of the form 

\^it)) = E Cn\un{t)) exp{-ient/h) (2) 



with time-independent coefficients c„. Therefore, under conditions of perfectly 
coherent evolution the occupation probabilities |c„p of the Floquet states 
remain constant, despite the periodic time dependence of the Hamiltonian. 
These coefficients then keep a memory of how the system has been prepared, 
that is, how the time-periodic forcing has been turned on in the past, whereas 
the Floquet states and their quasienergies are independent of this history. 

If, however, the evolution is not coherent, but the system is in contact with 
some sort of thermal environment or heat bath, what would be the equilib- 
rium distribution of the Floquet-state occupation probabilities which estab- 
lishes itself in the long run, irrespective of the particular initial condition? 
An important st e p tow ards a general answer to this question has been taken 
by iBreuer et al.l (120001 ) . who have studied anharmonic oscillators subjected 
to strong time-periodic forcing while being weakly coupled to thermal de- 
grees of freedom provided by a surrounding, and who have shown that the 
corresponding quasistationary density matrix actually is diagonal in the Flo- 
qu et representation. This line of inv estigation has recently been taken up 
by iKetzmerick and WustmannI (j2010l ). with a detailed view on the classical- 



quantum correspondence, after general asp ects of such "periodic thermody- 



namics" had been discussed by iKohnI (J200ll ). 



Turning now from the single-particle dynamics to that of an ideal gas of identi- 
cal quantum particles subjected to some periodically time-dependent external 
influence, such as a dilute gas of Bose particles stored in a confining poten- 
tial which is modulated periodically in time, an even more intriguing question 
poses itself: What are the statistics of the Floquet-state occupation numbers? 
In other words, what is the replacement for the familiar Bose-Einstein (or 
Fermi-Dirac, when dealing with Fermions) distribution function in a periodi- 
cally time-dependent setting? 

In order to app reciate what lie s behind this question, let us briefly retrace the 



steps which led lEinsteinl (119251 ) to what is now known as Bose-Einste in statis- 



tics, by now having become basic knowledge of statistical physics ( iPathria 



19961 ): Consider an isolated ideal Bose gas consisting of A^ ^ 1 particles pos- 
sessing a total energy -Etot, being kept in some large confinement. In view of 
the corresponding denseness of the single-particle spectrum we divide the en- 
ergy axis into small intervals, the ith such interval containing gi single-particle 
states with practically identical energy Ei. We then denote the number of par- 
ticles occupying the Qi states in the ith interval as n^, and ask for the most 
probable set {n*} of these cell occupation numbers. All admissible sets {rii} 
obviously have to comply with the two constraints 

J2n, = N (3) 



and 

J2niEi = Etot; (4) 



the latter holding approximately, reflecting the above coarse-graining proce- 
dure, but to good accuracy. The number of microstates associated with a given 
set {ui} is 



because the binomial coefficients quantify the number of possibilities to dis- 
tribute Hi indistinguishable particles over g^ states "with repetition", as is 
characteristic for Bosons. Assuming that both n^ ^ 1 and g^ ^ 1, and em- 



ploying Stirling's formula, one then obtains the entropy 



\nn[{n,}] = Y: 



n,ln(l + ^)+(7.1n - + 1 



(6) 



The searched- for most probable occupation numbers {^*} maximize this ex- 
pression, subject to the two constraints ([3]) and (jl]). Introducing a Lagrangian 
multiplier a to account for the conservation of the particle number, and a 
further multiplier /3 to account for the conservation of energy, we are thus led 
to the variational problem 



\ i i / 1 







(7) 



After brief calculation, this gives 



E 



\n[l + — ]-a-(3Ei 



Srij = 



(8) 



and thus results in the Bose-Einstein distribution 
1 



n, 



Qi exp{/3Ei + a) -I 



(9) 



quantifying the most probable occupation number of an individual single- 
particle state with energy Ei under microcanonical conditions. Finally, the 
multiplier /3 is identified with l/lk^T), with ^b denoting Boltzmann's constant 
and T being the temp erature of the gas, while a = —u/(k-p,T) is related to its 



chemical potential /x (JEinsteinl . 1 19251 : iPathrial . 119961 ) . 



When trying to adapt this reasoning to time-periodically driven quantum 
gases, and to determine the corresponding Floquet-state occupation numbers, 
we can immediately carry over the first constraint ([3]): The sum of all occu- 
pation numbers has to equal the total number of particles. This may seem 
trivial, but it is not, because already here we are making essential use of 
Floquet theory: It is only because the coefficients c„ in the expansion ([2]) 
are time-independent that we can assign occupation numbers to the Floquet 
states; only then can one ask for the associated distribution function. 

It is, however, not obvious what becomes of the second constraint (jl]), because 
energy eigenvalues do not exist for a time-periodic system, and the quasiener- 
gies of the Floquet states cannot be considered as their proper substitutes 
in this context. To understand the quasienergy concept more deeply, we in- 



sert the Floquet states ([T]) into the time- dependent Schrodinger equation, and 
deduce 

(H^'\t) - m^] \Un{t))) = BnKit))) . (10) 



The conspicuous notation employed here, i.e., writing the Floquet functions 
\un(t))) with a double right angle instead of \un(t)) as in Eq. ([T]), has a deep 
significance. Namely, Eq. ( ITOj) has to be regarded as an eigenvalue equation 
for the quasienergies; within the Floquet framework, this eigenvalue equa- 
tion takes over the role played for time-independent systems by the station- 
ary Schrodinger equation. Unlike the latter, this quasienergy eigenval ue equa- 



tion ( jTOj) lives in an extended Hilhert space of T-periodic functions (ISambe 



19731 ). In contrast to the usual physical setting, in which "time" emerges as an 
evolution variable, in that extended space time plays the role of a coordinate, 
and therefore needs to be integrated over in the associated scalar product: 
Denoting the momentary scalar product of two T-periodic functions \ui{t)) 
and \u2{t)) in the phy sical space as {ui{t)\u2{t)) , their scalar product in the 



extended space reads (jSambd . Il973l ) 



1 ^ 
{u,\u,)) ^ - j dt {u,{t)\u2{t)) . (11) 





Thus, we write \un{t)) when considering a Floquet function in the physical 
Hilbert space, as in Eqs. ([1]) and ([2]), whereas \un{t))) refers to that same 
function when viewed as an element of the extended space, as in Eq. ( fTOl) . 

This apparently formal observation has an important, physically meaningful 
consequence. Suppose that \unit))) is a solution to the eigenvalue equation (ITOl) 
with quasienergy e„, define u = 2tt/T, and let m be an arbitrary (positive, 
zero, or negative) integer. Then \un{t) exp{imut))) also is a T-periodic solution, 
with quasienergy £„ + mtiuj. For m 7^ these two Floquet functions obviously 
are orthogonal with respect to the scalar product flTTl) . and represent different 
eigenf unctions in the extended space. On the other hand, when returning to 
the physical space and forming the actual Floquet states ([T]), one has 

\un{t) exp(imco't)) expf — {{sn + mhu)t/h] = \un{t)) exp(— i£„t//i) , (12) 

so that both solutions represent the same physical state. Hence, a Floquet state 
does not correspond to a single solution to Eq. (fTOj) . but rather to a whole class 
of such solutions, labeled by the index n, while individual representatives of 
such a class are distinguished by the integer m. This subtlety already reflects 
itself in the expansion ([2]): Although all functions {\un(t) exp{imu)t)))} , with 



n ranging over all state labels and m extending over all integers, are required 
for the completeness relation in the extended space, only one representative 
from each class is required in Eq. ([2]), where no sum over the "photon index" 
m appears. By the same token, the quasienergy of a physical Floquet state 
is determined only up to an integer multiple of the "photon" energy tiw. In 
accordance with the analogous terminology used in solid-state physics, the 
quasienergy spectrum is said to consist of an infinite set of identical "Brillouin 
zones" of width hw, covering the entire energy axis, each state placing one of its 
quasienergy representatives in each zone. This means, in particular, that there 
is no natural "quasienergy ordering": Without additional specification it is 
meaningless to ask whether one given Floquet state lies "below" another. And 
there is still a further complication of a more mathematical nature: Because 
one finds one quasienergy representative of each state in each Brillouin zone, 
the quasienergy spectrum is "dense", and it may be technically difficult to 
decide whether one has a dense pure point spectrum, so that the expansion ([2]) 
can be taken literally and the system is stable, possessing a quasiperiodic 
wave function, or whether the re is an absolutely c ontinuous spectrum, allo wing 



for diffusive energy growth (JBunimovich et al.l . Il991t iHowlandl . Il992al ). We 



bypass this problem by restricting ourselves to stable systems with a pure 
point spectrum, such as th e forced anharmo nic oscillators with superquadratic 



potentials investigated by iHowlandl ( 1l992bl ). 



These considerations clearly reveal that, although there are Floquet-state oc- 
cupation numbers, the required analog of the energy constraint ^ for time- 
periodically forced microcanonical close-to-ideal Bose gases cannot involve 
uasienergies directly. H owever, each Floquet state possesses a mean energy 



Fainshtein et al.l . Il978[ ) 



1 '^ 

En = -jdt{Un{t)\H^^\t)Wn{t)) = ((«„ |iJ« l^n)) (13) 





which obviously is independent of the choice of the representative !«„)). If we 
now repeat the coarse-graining procedure employed in the formulation of the 
energy constraint (jlj), that is, if we divide the energy axis into small cells, the 
ith cell now containing gi Floquet states having almost the same mean energy 
Ei, and if {rij} denotes the corresponding sets of cell occupation numbers, 
then it seems reasonable to demand 

^nji^j = const. (14) 



for an isolated system, thus obtaining a constraint which effectively may re- 
place the previous Eq. (j4]). Here we assume that each driven particle sees the 
other driven particles as its heat bath, and undergoes subsequent short relax- 



ation events only on a time scale which is large compared to the period T. 
In general, the justification of this constraint may require the consideration 
of particular physical setups, and some more specifications may be needed. 
But if we tentatively accept this constraint f ll4p as it stands, it is immediately 
clear how to proceed: We assign the entropy ([6]) to a time-periodically forced 
isolated Bose gas with Floquet-state occupation numbers {rii}, account for the 
particle number conservation (|3]) again by introducing a Lagrangian multiplier 
a, and incorporate the conservation of the mean energy flT4|) with the help of a 
further Lagrangian multiplier 7. Then the most probable set of Floquet-state 
occupation numbers {n*} is determined by the variational equation 

s(\nn[{ni}]-aY.n,--fY,n,^A =0. (15) 



Be cause thi s prob lem exactly parallels the problem (^^ already considered 



by lEinsteiru (119251 ). we can directly take over the solution: Under the above 



conditions, the most probable set of Floquet-state occupation numbers is 

n* 1 

Qi expi^jEi + a) - 1 ' 

suggesting that 

f(En; a, 7) = . ^^^ , , (17) 

exp(7^„ + a) - 1 

is the expected occupation number of an individual Floquet state with mean 
energy En- Of course, this reasoning can easily b e adpated to Fe rmions, leading 



to the familiar "plus" sign in the denominator (jPathrial . Il996l ) . 



These deliberations may have some profound physical implications. To begin 
with, a time-periodically forced, stable Bose gas manifestly does not con- 
stitute a nonequilibrium system, as one might naively assume, but rather an 
equilibrium one, with equilibrium parameters a and 7. Most noteworthy, these 
parameters depend on the form of the periodic time- dependence, because so 
do the Floquet states, and hence their mean energies flT3|l . which enter into 
the modified constraint flT^ . Therefore, when writing 



fcRfc) 



in analogy to the familiar relation /3 = l/lk^T), the parameter B introduced 
here is a temperature-like quantity which can be varied by changing, e.g., 
the strength of an external time-periodic force; in general, B will differ from 



the usual temperature T the gas would have if there were no such forcing. 
Likewise, when setting 



z/ is the corresponding, forcing- dependent chemical potential. 

It follows that there are Floquet condensates: If the 0-temperature is suffi- 
ciently low, a time-periodically forced close-to-ideal Bose gas condenses into 
the single-particle Floquet state possessing the lowest mean energy. This de- 
duction immediately suggests a further possibility: Suppose that the Hamil- 
tonian under consideration has the natural form H^^\t) = Hq + Hi(t), where 
Hi{t) represents a time-periodic force acting on the unperturbed system Hq 
with adjustable strength. In the regime of perturbatively weak forcing the 
mean-energy ordering of the Floquet states of H^^^ (t) is likely to be the same 
as the energy ordering of the eigenstates of Hq from which they have devel- 
oped, but in the nonperturbative regime this is no longer guaranteed; here the 
Floquet state possessing the lowest mean energy not necessarily is connected 
to the Ho-ground state. Under such conditions the Floquet state into which 
an ultracold, time-periodically forced Bose gas condenses is determined by the 
strength of the forcing, and may change when the latter is varied. 

We close this section by emphasizing that the above reasoning hinges on two 
key issues: The assignment of occupation numbers to the Floquet states, which 
is borne out by theory, and the adoption of the mean-energy constraint fll4p as 
a building principle for a microcanonical Floquet ensemble, which may require 
further thoughts. Important pieces of evidence supporting our theoretical de- 
ductions are provided by the experimental observations reported in Sec. [6] 
below. In order to facilitate the interpretation of our measurements, we first 
lay out some important ingredients in Sees. [3] to O 



3 The experimental setup: Shaken optical lattices 



A one-dimensional (ID) optical lattice potential, which we write in the form 

Vn 

V{x) = — cos(2A;lx) , (20) 



is generated by two linearly polarized counterpropagating beams of laser ra- 
diation w ith wave number fcr,, its dept h Vq being proportional to the laser 



intensity (JMorsch and Oberthalerl . l2006l ). Here the lattice coordinate x refers 



to the laboratory frame of reference. Considering an atom of mass M moving 



10 



in this potential, the lattice depth Vq is measured in units of its single-photon 
recoil energy E,., 



Er 



~2M 



(21) 



typical optical lattices are up to tens of recoil energies deep. For example, 
when working with ^^Rb in an optical lattice made up from laser light with 
wavelength A = 27c/ki^ = 852 nm one obtains E,- = 1.31 ■ 10^^^ eV. Thus, the 
typical depth of an optical lattice is 11 orders of magnitude lower than that 
of the lattices encountered in traditional solid-state physics. This implies that 
the characteristic frequencies of ultracold atoms in optical lattices fall into the 
lower kilohertz regime. 



The pairs of counterpropagating beams needed to create optical lattices in 
one, two, or three dimensions can either be realized by splitting a laser beam 
in two, or by retro-reflecting a beam off a mirror. In the former case, a small 
frequency difference Az/(t) introduced between the two splitted laser beams 
with the help of aco usto-optic r nodu l ators makes the lat t ice ni o ye with velocity 



v(t) = Au(t)X/2 fiNiu et al.l. Il996 : iBen Dahan et all . Il996l : Madison et al 



19981 : iLignier et al.l . 120071 : ISias et al.l . 120081 ) . Therefore, in the laboratory frame 
of reference an atom then experiences the potential 



Vii,hix,t) 



YCos(2A;l[. 



X 



Mt)] 



(22) 



where 



Xoit) = ^JdTAu{r). 



(23) 



When using retro-refiected beams, lattice motion can be achieved by mounting 
the mirror on a piezo-elec tric actuator, so tha t it shifts acco rding to some 
prescribed protocol Xo(t) (IZenesini et al.l . l2009l : lAlberti et al.l . |2P09). As the 
shift range of such devices is in the micrometer regime, this method is suitable 
for inducing an oscillatory motion as needed for shaking lattices, but not for 
applying a constant acceleration. 



In the following steps, the Schrodinger equation for the translational motion 
of an atom in the shifted lattice, governed by the Hamiltonian 



H\a.h{t) 



2M 



+ V[ab(a;,t) 



(24) 



11 



is subjected to a unitary transformation to the frame of reference co-moving 
with the lattice. First, the corresponding shift in position is implemented by 
means of the unitary operator 

t/i=expQxo(t)p) , (25) 



implying 

UixUl = x + Xo{t) (26) 



and 



^'(-4)^* = -4-^»W''- 



(27) 



thus leading to the new Hamiltonian 

^(^) = 9Xf (^ " MXo(t))' + ^ cos(2A:lx) - ^Mt)' ■ (2^ 



Next, the shift of momentum showing up here is compensated through 

U2 = exp(^-^MXo{t)x^ , (29) 



givmg 

U2pUl=p + MXo{t) (30) 

together with 

U2 (-i^^) U^2 = -i^^ + MXo(t) X . (31) 



Finally, the merely time-dependent energy shift appearing in H{t) is removed 
by the operator 

^3 = expU^|^|drXo^(r)j , (32) 



12 



since 



^4-4)''--4+f-*»w- p^' 



The full transformation to the co-moving frame is then realized by the com- 
bined operation U = USU2U1, giving 

[/(i7iab(t) - ihdt)u^ = Ho + H^{t) - ihdt , (34) 



with 

Ho = ^ + ^cosi2k^x) (35) 



denoting the single-particle Hamiltonian pertaining to the undriven lattice, 
and 

Hi{t) = -F{t)x (36) 



introducing a homogeneous inertial force F{t) acting in the co-moving frame, 
determined by the lattice motion according to 

F{t) = -MXo{t) . (37) 



In particular, when a purely sinusoidal frequency shift with an initial phase 
is suddenly applied at t = 0, 

,0 , t<0 , 

Au{t) = { (38) 

Az/max sin(wt + 0) , t > , 



Eqs. ([22]) and (|23D give 

Mab(a;, t) = Y cos(2A;L[a; + L{cos{ut + 0) - cos(0)}]) (39) 



for t > 0, so that the lattice is shaken in the laboratory frame with the 
amplitude 

L = ^^^^ . (40) 



13 



According to Eq. (123|) . one has Xo(t) = AAz/(t)/2. Therefore, this protocol ( 138|) 
effectuates a sudden jump of the lattice's velocity at time t = 0, unless = or 
(j) = 71. By virtue of Eq. fl37p . the force felt by the atoms in the co- moving frame 
of reference then is composed of a monochromatic oscillating drive acting for 
t > 0, and of a delta-like kick acting at the moment of turn-on: Formally, one 
obtains 

F(t) = -Fi cos(a;t + 0) Q{t) - — sin(0) 6{t) , (41) 

where 9(t) denotes the Heaviside function, and the driving amplitude is given 
by 

Fi = M^Az/^axW = MLoj^ . (42) 



Needless to say, in reality the delta-kick experienced at t = has a finite 
sharpness, determined by the short-time details of the actual velocity jump. 
Nonetheless, this kick has an important experimental consequence, as will be 
discussed later. 

As a dimensonless measure of the shaking or driving strength we introduce 
the quantity 

F,d _7T^U L 

nio 2 ujj- a 



where d = A/2 is the lattice constant. Thus, with driving frequencies on the 
order of the recoil frequency Uj. = Ej-fh, and modulation amplitudes (HOj) on 
the order of the lattice constant, one can easily reach the nonperturbative 
regime Kq > 1. It would be quite hard to realize corresponding conditions 
in laser- irradiated traditional solids without introducing, e.g., additional po- 
larization effects, or even damaging the sample. Therefore, strongly shaken 
optical lattices may also be viewed as "strong-field simulators" which al- 
low one to study even superstrong-field-induced multiphoton-like processes 
i n periodic potentials, su c h as interband transitions, in their purest form 



(JArlinghaus and Holthausl . 120101 ). 



4 The driven Bose-Hubbard model 



In principle, the time-periodic motion (139|) always induces transitions between 
the unperturbed energy bands of the optical lattice. However, if the energy 
scale huj associated with the modulation frequency remains small compared 
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Fig. 1. Exact energy gap A between the lowest two bands of an optical cosine 
lattice (j2Up . measured in multiples of the recoil energy Et- (full line), in compari- 
son with the shallow-lattice estimate (|44p (long dashes), and with the deep-lattice 
estimate (US} (short dashes). 

to the energy gap A between the lowest two bands, and if the driving ampli- 
tude remains sufficiently low, the dynamics of driven ultracold atoms remain 
restricted to the lowest band at least to good approximation. Since the band 
structure of a cosine latt ice ([20]) is d etermined by the characteristic values 
of the Mathieu equation (ISlaterl . Il952l ). the known expansions of these values 
( lAbramowitz and Stegunl . Il965l ) can be employed for obtaining estimates of 
the gap width, resulting in 



A/E, ^ Vo/i2E,) 



(44) 



for shallow lattices, and 



A/E, ^ 2JV0/E, - 1 



(45) 



for fairly deep ones. Fig. [T] demonstrates that these two approximations indeed 
provide a reasonable estimate of the exact band gap for all Vo/Ej., if one 
switches from the shallow-lattice result (jUj) to the deep-lattice formula (H5l) 
at Vq/E^ ^ 11.7. As a figure of merit, the exact gap width is A = 4.572 E^. for 
a lattice with depth Vq = lOi^r- 



In our experinaents (ILignier et al.l . 120071 : ISias et al.l . l2008l : lEckardt et al.l . 12009 



Zenesini et al.l . |2009| ) we work with Bose-Einstein condensates of ^^Rb con- 
sisting of about 5 X 10'^ atoms in a shaken ID optical lattice (A = 852 nm or 
842 nm) with depths Vq/E^ ranging from 3 to 10, employing scaled shaking 
frequencies hoj/Ej. between roughly 0.1 and 2 (corresponding to u/{2'k) be- 
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tween 0.3 and 6 kHz). Under such conditions the single-band approximation 
is viable, provided the forcing is not too strong, si nce crossing the gap would 
requi re higher-order multiphoton-like transitions ( lArlinghaus and Holthausl . 
2OIOI ). In order to model these driven, interacting many-body systems we re- 
sort to the cq - moving frame of refererice, arid then follow the s tandard route 



( jjaksch et al.l . Il998l : iJaksch and ZoUerl . l2005l : iBloch et al.l . |2008| ) . Employing a 



basis of site-localized Wannier functions pertaining to the lowest Bloch band, 
and assuming a reasonable lattice depth Vq/ E.^ ^ 1, it suffices to retain only 
the hopping matrix element J connecting neighborin g sites. The accura cy of 
thi s approximation has been assessed quantitatively by lBoers et al.l (120071 ) and 
by lEckardt et al.l (120091 ): for instance, when Vq/E^ = 10 the neglected matrix 
element connecting next-to-nearest neighbors is actually smaller than J by 
a factor of about 0.012. Moreover, the van d er Waals length of alkali atoms 
typically amounts to just a few nanometers ( JBloch et al.l . l2008l ) and thus is 
significantly smaller than the lattice constant d = A/2, so that only the on-site 
interaction among the atoms has to be accounted for; this is done in terms of 
a parameter U which quantifi es the interaction energ y of one pair of atoms oc- 



cupying the same lattice site ( ISchneider et al 

these three approximations (single-band, nearest-neighbor, and on-site) 



2OO9I ). In short, when adopting 

an 



ultracold gas of bosonic atoms in a time-periodically shaken optical lattice is 
described in th e co-moving fra r ne by the explicitly time-dependent rn any-body 
HamiltoniaiH dEckardt et al.l . l2005bl : ICreffield and Monteirol . I2OO6I ) 



Hit) = Ho + H,it) 



(46) 



where 






(47) 



is the standard Bqse-Hubbard model for a ID lattice (iFisher et al.l . Il989l : 
Jaksch et al.l . Il998l : IJaksch and ZoUerl . |2005| ). with b^ denoting the bosonic 



annihilation operator for atoms occupying the Wannier state at the £th site; 



ni 



h\h(^ is the number operator for that site. Recalling that the position 
operator x mediating the inertial force built into the single-particle Hamilto- 
nian (l36|l translates, on the many-body level, into d^i^h^, with d = A/2 for 
the distance between two sites, the time-periodic forcing now is introduced 
through 



Hi{t) = Kcos{ujt)J2^ne 



(4J 



We use the "hat" -symbol to indicate operators acting on the bosonic Fock space 
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where the strength K is related to the forcing amphtude ( H2|) hy K = Fid. 
Here we disregard an initial phase 0, which would refer to a sudden turn-on 
of the force, as in Eq. (1381) . This driven ID Bose-Hubbard rn odel (H^ may be 
extended to higher spatial dimensions ( IZenesini et al.l . |2009| ) ; we remark that 
approximate expressioi is relatin g the p aramete rs J and U to the l attice depth 
have been provided by IZwergerl (120031 ) and by iBloch et al. 



mm . 



5 Interference patterns produced by Floquet states 



It is essential to observe that the existence of Floquet states ([T]) hinges solely 
on the periodicity of the given Hamiltonian in time. Thus, besides the single- 
particle Floquet states which form the basis of the statistical considerations 
in Sec. El there also are Floquet states for periodically time-dependent, in- 
teracting many-body systems, such as periodically driven Rose gases, which 
incorporate both the periodic time-dependence and all interaction effects. 



5. 1 Signatures of interacting shaken Bose gases 



Here we are concerned, in particular, with the Floquet states of the driven 
Bose-Hubbard model (B^ . and their experimental signatures. In order to solve 
the quasienergy eigenvalue problem (fTOj) for this model, the quasienergy op- 
erator H{t) — i hdt h as to be diagonalized in the associated extended Hilbert 
space (ISambd . Il973l ). as reviewed in Sec. El To this end, let {n^} denote an 
admissible set of site-occupation numbers n^. The physical many-body Hilbert 
space then is spanned by the set of all Fock states 



iw)=n4SKac) 



'ni>\ 



(49) 



where |vac) is the "empty-lattice" state. A possible basis of the extended space 
would be provided by the products |{7^^}) exp(i?7iwt) with integer m. However, 
it is more useful here to employ the basis of Floquet-Fock states (or "dressed 
Fock states") given by 



\{ni},m)) 



Ijn^}) exp< —It— sin(c<;t) ^£n^ + imut > 
\ hw ^ J 



(50) 



which already diagonalize Hit) —ihdt for vanishing i nterw ell hopping strength 
J = (JEckardt et al.l . l2005bt : lEckardt and Holthaud . 120071 1 . Invoking the scalar 
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product (ITT!) , in this basis the matrix of the quasienergy operator has the 
structure 



{{{n'f,}, m'\ H{t) - ihdt |{raj, m)) = 5.m',m {W(}\ H^s + mhu |{nj) 

+ a-S^m',m){{n',}\V\{ne}), (51) 

where Hes is a time-independent Bose-Hubbard Hamiltonian of the famihar 
form (H7I) . but with the hopping matrix element J multiphed by a zero-order 
Bessel function Jq with the argument Kq = K/{hu), resulting in the effective 
hopping strength 



Jeff = JJoW(M 



(52) 



The operator V, which is off-diagonal with respect to the photon index m, 
contains nearest-n eighbor couplings J scaled by higher-order Bessel fu nctions 



±Jm-m{K/{huj)) flEckardt et al.l . l2005bl : lEckardt and Holthausl . 120071 ). Thus 



the diagonalization of the quasienergy matrix 05 ip constitutes a tremendous 
problem of a kind not usually considered in traditional many-body physics: 
Infinitely many ffgff-blocks, each one corresponding to the full Hamiltonian 
matrix of a time-independent Bose-Hubbard model (H7j) with modified hop- 
ping matrix elements (152|) . are shifted against each other in energy by integer 
multiples of hu, and are coupled by V"-blocks. The latter blocks, in their turn, 
embody hopping elements which are multiplied by Bessel-function factors with 
indices reflecting the distance of the respective block from the main diag- 
onal. In effect, these V^-couplings cause multiphoton -like resonances among 
the st ates described by the shifted diagonal blocks (JEckardt and Holthausl . 
2008bl ). While the physics of this problem has not yet been explored in full 
generality, its high-frequency regime is comparatively transparent. Namely, if 
hw is large compared to the two energy scales J and U (while remaining still 
small compared to the band gap A of the underlying optical lattice, so that 
the single-band treatment remains applicable), it is a good approximation 
to neglect the couplings induced by V altogether, so that the driven sys- 
tem (H^ reduces to an undriven system described by Hsq, meaning that the 
effect of the time-periodic force essentially is to "renormalize" the hopping 
matrix element according to Eq. fl52|) . With the V^-couplings out of action, 
all ifeff-blocks are equivalent, and returning from the extended to the physi- 
cal Hilbert space is tantamount to considering only one such block from the 
outset. Since the ratio J/U governs the s uperfiuid-to-Mott i r isulator t ransi- 



tion intrinsi c to th e Bose-Hubbard model (iFisher et al.l . Il989l : IZwergerl . 12003 



Bloch et al.l . l2008l ). the renormalization of J in response to high-frequency 
forcing implies that it is possible to induce that transition through adjusting 
the paramet ers of the driving force , while keeping the lattice dept h and hence 



U constant (JEckardt et al.l . l2005bl : lEckardt and Holthausl . |2008a|). This type 



of coherent control over the superfluid-to-M ott insulator t r ansiti on has been 
demonstrated in a pioneering experiment by lZenesini et al.l (120091 ) . It involves 
adiabatic following of the i^eff-ground state when the driving amplitude is 
slowly changed; this has meanwhile been studied i n detail with the help of 
numerical simulations by iPoletti and KoUathl ( 1201 ll ). 



As in the time-independent case, experimental information about the many- 
body state \ip{t)) is obtained by time-of-flight absorption imaging: Switching 
off the time-periodically shaken lattice potential at some moment t, then let- 
ting the matter wave expand (neglecting interaction effects during that ex- 
pansion), and finally recording its density in spac e after a sufficiently long ex- 



pansi on time yields the momentum distribution (IZwergerl . l2003l : iBloch et al. 
2008h 



n{p,t) = {ij{t)\d\p)a{p)m)) 



\w 



(p)P^exp Q(r - s)pci) (^(t)|6t6J^(t)) 



(53) 



where a{p) is the annihilation operator for a free-particle state with momen- 
tum p in the direction of the lattice, and w{p) is the Fourier transform of the 
Wannier function pertaining to the lowest Bloch band; again, d = A/2 denotes 
the lattice constant. Thus, apart from the factor |w(p)p the observed momen- 
tum distribution is given by the Fourier transform of the system's one-particle 
density matrix {'ip{t)\blb^\ip{t)). Let us now assume that for a given driving 
amphtude K/{hu), and hence for a given value of Jes, an eigenstate of Hes 
with energy e takes the form 



li'esit)) =exp( 



-iet/h) ^7{nalW}) 

{ne} 



(54) 



with certain coefficients 7{„j. Under conditions such that the high-frequency 
approximation detailed above is valid, the corresponding many-body state in 
the driven optical lattice is then obtained by replacing the Fock states (I49p in 
this expansion (1M|) by the Floquet-Fock states (ISUl) with m = 0, say, resulting 



m 



\m) 



exp{-iet/h) Y, 1{n,}\{ni}) expl -i— sin(wt) ^£n^ 

{rti} ( £ 



(55) 



In order to evaluate the momentum distribution ( 153|) for this particular many- 
body Floquet state (l55l) we then compute 
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{m\blbs\m)= E lw,}^{nMWe}\blbs\{ne}) 

{ne},{n'} 



f K 1 

exp< — i^" sin(ci;t) ^ i{ne — n'^) > . 
[ nw f, J 



X exp\ —i^r- sin(ci;t) y^ i{ne — n'f) } . (56) 



Now the only nonvanishing matrix elements are those with rii = n'^ for £ ^ r,s] 
together with nj, = n.^ + 1 and n'^ = Ug — 1. This implies 

J2i(ne-n',) = -ir~s), (57) 

e 

which allows us to take the phase factor appearing on the right-hand side of 
Eq. (!56|) out of the sum, giving 

imiKK \m) = (^effl bib, \^Pes) exp|i(r - s)^ sin(u;t)| . (58) 

Inserting this result into the representation flS5]) . we immediately obtain the 
momentum distribution provided by a single Floquet state of the driven Bose- 
Hubbard model (H6|) in the high-frequency regime: 



n(p, t) = |w(p)p E(^eff I bib, IV^eff) 

X exp< — (r — s) ( pd H sin(a;t) | > , (59) 

where, by assumption, \ipes) is the associated energy eigenstate of ^eff- 

However, an important step is still missing in order to connect theory with 
experiment: The driven Bose- Hubbard model, with the inertial force being 
incorporated through the driving term (B8|) . refers to the co- moving frame of 
reference as considered in Sec. [3l whereas measurements usually are performed 
in the laboratory frame. Thus, in order to obtain the momentum distribution 
niab(piab5^) as recorded by an observer in the laboratory frame, we still have 
to invert the transformation (I5U]) . This is done in the general case by writing 

p = piab - MXo(t) , (60) 

and through 

P = Pish -1 sm(a;t) (61) 

for the driven Bose-Hubbard model with forcing (1481) . as corresponding to a 
frequency variation Ai^(t) = Az/inaxSin(a;t) between the counterpropagating 
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laser beams, or to lattice motion Xo{t) = Lujsm{ujt), keeping in mind the 
relation (H2|) for the driving amplitude, together with the definition K = Fid. 
Therefore, we finally have 



^lab(Plab, t) = |w(piab " K / {ud) Sm{ujt) 



X 



^exp f -(r - s)piab'ij (V^cffl hlh^ |V^eff) • 



(62) 



This is quite a significant observation: The momentum distribution of a mat- 
ter wave occupying a many-body Floquet state in a time-periodically shaken 
optical lattice ( 139|) in the high-frequency regime equals that of the associated 
energy eigenstate of -ffeff? which is obtained from the undriven Bose- Hubbard 
Hamiltonian f H7|) through replacing the hopping matrix element J by Jeg 
as defined in Eq. flS2]) . insofar as exactly the same Fourier transform of the 
one-particle density matrix appears in both cases. Hence, even though the 
position of the lattice is periodically shifted, the peak pattern observed in the 
laboratory frame does not move. The effect of the time-periodic shift is seen 
only in the envelope of that pattern, given by the Fourier transform w of the 
Wannier function, the argument of which is modulated periodically in time in 
accordance with Eq. f l6T]) . Thus, apart from this modulation the experimental 
signature of the superfluid-to-Mott insulator transiti on occurring in shak en 



optical lattices upon changing the driving amplitude (jZenesini et al.l . |2009| ) is 



the same as that of the tra nsition o c curing in a stationary la ttice in response 



to a variation of its depth (jZwergerl . l2003l : iBloch et al.l . l2008l ) 



It may be useful to point out that the experimental signatures differ from 
the above description if the optical lattice is not shaken, but kept at rest, 
while the force is induced by means of the time-periodic modulation of a 
lev itation gradi e nt wh ich "stirs" the condensate, as done in the experiments 
by iHaller et al.l (120 lOl ). In this latter situation the model ( HHj) actually de- 
scribes the driving force in the very reference frame in which the momentum 
distribution is recorded, so that one obtains an oscillating interference pat- 
tern given di r ectly b y Eg. (I59l). In con t rast, in the ex p erime n ts reported by 



Lignier et al.l (120071 ). ISias et al.l (120081 ). lEckardt et al.l (120091 ) . IZenesini et al 



( 120091 ) . and in Sec. [6] below, the lattice is shaken in the laboratory frame 
according to Eq. ( 151?]) . Hence, here the "micromotion" is taken out of the in- 
terference pattern by means of Eq. (l^Tl) . which connects the momentum p in 
the co-moving frame to the momentum piab observed in the laboratory frame. 
The fact that the resulting interference pattern does not move in the labo- 
ratory frame (JEckardt et al.l . l2005bl : lEckardt and Holthausj . l2007l . l2008al ) may 
sometimes facilitate its interpretation. 
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5.2 Signatures of ideal shaken Bose gases 



With respect to the Floquet condensates envisioned in Sec. [21 the case of an 
ideal Bose-Einstein condensate in a shaken optical lattice now is of particular 
interest. For vanishing interaction U = the time- dependent Schrodinger 
equation for the model f H6|) can easily be solved exactly for any type of forcing 
described by a driving term 

H,{t) = -F{t)dJ2ehe, (63) 



without requiring a specific time-dependence of the force. Therefore, we aban- 
don time-periodic forces and the Floquet picture for the moment being, and 
consider an initial-value problem instead: We assume that F{t) vanishes for 
t < 0, and is switched on sharply at t = 0, but is arbitrary otherwise. Gener- 
alizing the previous Eq. (1411) . we thus impose a force 

F{t) = -MXo(t) 0(t) - MXo(0+) S{t) , (64) 



where the second term accounts for the sudden velocity jump of the lattice 
from Xo(0— ) = to an arbitrary value Xo(O-l-) at t = 0. Moreover, we con- 
sider a lattice with M^ ^ 1 sites and disregard finite-size effects, so that the 
operator 

4(0) = ^Eexp(i^^ci)S{ (65) 



creates a particle in the Bloch state with quasimomentum hk. We then assume 
that the initial state at t = is an ideal A^-particle condensate occupying such 
a Bloch state. Even if A; = might be the only experimentally realistic option 
here, we do not impose this restriction at this point. For i > 0, after the 
force has been turned on, the resulting A^-particle wave function can then be 
written in the form 

1 r + iN 
\Mt)) = -T^Jclit)] |vac). (66) 



Here the creation operator c|,(t), given by 

cl(t) =exp[-^|drE(gfc(r)) I ^ ^exp(i%(t)d) SJ , (67) 
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refers to a so-called Houston state, also known as accelerated Blocli state 



( jHoustonl . Il940l : lEckardt et al.l . 120091 ) . This nomenclature stems from the fact 
that the time- dependent wave number qk{t) appearing here has to obey the 
"semiclassical" acceleration law 

hqk{t) = F{t) . (68) 



Therefore, using the particular connection (1641) between the inertial force and 
the lattice motion, and requiring that qk{t) be equal to the wave number k of 
the initial state for t < 0, we have 



1 } 
qk{t) = k + -JdTF{T) 





= k - y (Xo(t) - Xo(0+)) - -Xo(0+) 

= k-^Mt) (69) 

for t > 0, having properly accounted for the delta- kick at the moment of 
turn-on. Finally, the expression 

E{k) = -2Jcos{kd) (70) 



appearing in the exponential of Eq. (!67j) denotes the single-particle dispersion 
relation describing the energy band provided by the Hamiltonian (HTIl when 
U = 0. 

For such a noninteracting X-particle wave function fIBB]) the momentum dis- 
tribution (!53|) is given, apart from the factor |w(p)p, by 



^exp(i(r - s)pd/h) {Mt)\blk\Mt)) 

r,s 

■ NY,exp{ir[p/h-qk{t)]d). (71) 



The return to the laboratory frame now is achieved with the help of Eq. (160|) . 
and the expression (!69|) for qk{t): Using these, one obtains 

p/h - qk{t) = puh/h - k . (72) 



Therefore, the interference pattern provided by an ideal Houston condensate in 
time-of-fiight absorption imaging again is stationary in the laboratory frame. 
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and the peak positions are not affected by the force: Even after the force ([61 
has been switched on, the interference peaks are still permanently centered 
around the wave numbers k mod {27r/d), as they had been for the unforced 
initial condensate. It needs to be stressed that this peculiar feature is crucially 
dependent on the delta- kick which accompanies the sudden turn-on of the force 
in the co- moving frame. 

In particular, let us now consider a monochromatic oscillating force switched 
on instantaneously at t = with starting phase according to Eq. (HT1) . which 
leads to 

gfc(t) = fc - -— sin(wt + {/)) (73) 

lieu 



for t > 0. In this case there is a close relation between the A^-particle Houston 
states ( l66l) . which have been constructed for t > as solutions of an initial- 
value problem, and A^-particle Floquet states, which presuppose a perfectly 
time-periodic force F{t) = —Fi cos{ut + 0) acting at all times t. In order to 
obtain these Floquet states, one only has to extend qk(t), as given by Eq. (!73|) . 
to all t, and use this expression in Eq. ( l67l) . now to be considered for all t. 
Then the emerging "extended" Houston states (1661) have the basic form ([T]): 
By construction, qk(t) is T-periodic, with T = 27i/uj, but the energy integral 
in the exponent of Eq. fl671) is not. However, writing 



exp --ydrF(gfe(r)) 

= exp --^ Jdr [E{qk{T)) - e{k)] J exp ( - ie{k)t/h) 
with the help of the one-cycle-averaged energies 

T 



£(fc) = i|dri?(g,(r)) 





= -2J Joi ^] cos(kd) , (74) 

\nio J 

the first of these exponentials is T-periodic and therefore part of the Floquet 
functions, whereas the occurrence of the averages (17^ in the second expo- 
nential, accompanied by the time t, allows one to identify them as quasiener- 
gies. The Houston-Floquet states thus found are particularly simple examples 
of spatiotemporal Bloch waves, incorporating both the spatial periodicity of 
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the lattice and the temporal periodicity of the driving force on equal foot- 
ing. They are labeled by the same quantum numbers k as the customary, 
time-independent Bloch waves to which they reduce in the absence of the 
drive, while their time-evolution, apart from the time-periodic motion incor- 
porat ed into the moving wave number s qkit), is specified by the quasienergies 
6{k) (jArlinghaus and Holthaud . l201ll ). Observing that one again encounters 
here the effective hopping matrix element Jes introduced in Eq. fl52l) . the 
quasienergy-quasimomentum dispersion relation takes the form 



e{k) = — 2 Jcfj cos{kd) mod hw 



(75) 



which differs from the original dispersion relation (!70l) of the undriven lattice 
only through the replacement of J by Jck- That same replacement had been 
met before in the context of the interacting many-body system described by 
the driven Bose-Hubbard model (1461) . when constructing the approximate ef- 
fective Hamiltonian ffeff pertaining to the high-frequency regime. In contrast, 
in the noninteracting case considered here no approximations have been made; 
Eq. fl75|) holds exactly for all driving frequencies. 



The fact that the quasienergies fj75l) for the driven, noninteracting Bose- 
Hubbard model could be calculated by taking the time-averages (TMl) further 
suggests that these quasienergies do coincide with the corresponding mean en- 
ergies fITS]) . as considered in Sec. [2J Indeed, taking a single-particle Houston- 
Floquet state iV'fc(^)) = \ukit)) exp ( — ie{k)t/h) as constructed above, one 
easily confirms the identity 



E{k) = {{uk\H{t)\uk)) = -2J,scos{kd) 



(76) 



so that here the mean energies of the Floquet states actually equal the quasi- 
energies fl75l) . disregarding their "mod /io;" -multiplicity. In the case of a gen- 
eral T-periodic single-particle Hamiltonian H^^'^ (t) with Floquet states ([1]) this 
is not the case, since then 



En = {{Un\H^'^^ - lhdt\Un)) + ((n„|i/l<9, 
= £„+ {{Un\lhBt\Un)) , 



■t\Un 



(77) 



as discussed in detail by iFainshtein et al.l ( 1l978l ). 



Interestingly, the quasienergy band (175]) "collapses" w hen Kq = Fidjihw) 
equals a zero of the Jo Bessel function ( jHolthausl . Il992l ). implying that an ar- 
bitrary single-particle wave packet driven under such conditions cannot spread, 
but reprodu ces itself T-periodically. T his phenomenon, termed "dynamic lo- 
calization" (JDunlap and Kenkrd . 119861 ). has recently bee n observed w i th di- 



lute Bose- Einstein condensates in driven optical lattices (JLignier et al.l . 12007 
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Eckardt et al.l . |2009[ lArlinghaus et al.l . 1201 ll ). With a view towards future 
applications, the dependence of the quasienergy band width on the driving 
strength also has been identified as a me ans of controlling transport in sys- 
tems with attractive pairing interactions ( jKudo et al.l . |2009| ) . 



Once again, it is instructive to compare the above results, obtained for a shaken 
optical lattice, to the corresponding physics when the lattice is at rest, while 
the condensate is stirred by a harmonically modulated levitation gradient 
( jHaller et al.l . I2OIOI ). Then it is actually possible to impose an instantaneously 
turned-on oscillating force 



F{t) = -Fi cos{ujt + 0) 0(t) 



(7J 



without the additional delta- kick present in Eq. (I^Tl) . which necessarily ap- 
pears when an inertial force is abruptly turned on. The moving wave numbers 
Qkit), which previously had been determined in Eq. (169!) as solutions to the 
equation of motion (!68l) . now have to be replaced by the solutions q(t) to the 
corresponding equation hq(t) = F(t), giving 





Fi Fi 

= k + —- sin(0) — —- sinfwt + 0) 
niij njjj 

for t > 0, instead of Eq. fl75]) above. Therefore, one formally has 
q{t) = qk+Ak{t) 



(79) 



^0) 



with a wave-number shift 

Fi 
Ak = -—sinicj)) , 



(81) 



which means that an additional momentum hAk is imparted on the conden- 
sate particles when the stirring force is switched on rapidly, as described by 
Eq. (1781) . Thus, besides the micromotion the resulting time-of-flight absorption 
images would also show a shift of the peak positions. Both of these features 
are absent in experiments with shaken lattices: It really matters whether the 
condensate is "shaken" or "stirred" ! 



We note that iKudo and Monteirol (l2011al Jbl). focusing on stirring forces of the 
form (ITHl) . have suggested to introduce effective dispersion relations which 
also incorporate the initial phase 0. While one is free to adopt this viewpoint. 
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it seems to obscure the conceptual simplicity of the Floquet approach: The 
quasienergy band (175|) characterizes the time-periodically driven system as 
such, regardless of the way the drive has been switched on, that is, independent 
of the phase which parametrizes the sudden turn-on (178 p . or of any other 
parameters which specify other, equally possible turn-on protocols. Because 
this quasienergy band fl75|) actually consists of eigenvalues of the quasienergy 
operator it even provides a dispersion relation in the usual sense of solid-state 
physics which allows one to compute group velocities by taking its derivative, 
properly ev aluated at that quantum number k around which the wave packet 



is centered f lArlinghaus and Holthausl . 1201 ll ) 



With regard to shaken condensates, we may summarize our considerations 
as follows: If a perfectly ideal Bose-Einstein condensate in an optical lattice 
initially occupies a Bloch state with wave number k, and then is abruptly being 
shaken and thus subjected to the force (HT!) under single-band conditions, it 
permanently populates a single Floquet state labeled by the same wave number 
k, regardless of the phase 0. 



6 Experimental findings 



In order to substantiate the relevance of the above ideal-gas considerations 
for laboratory experiments, we took time-of-fiight absorption images of dilute 
®^Rb-condensates released from shaken ID optical lattices fl39l) . employing 
various driving frequencies and amplitudes. In Fig. [2] we display typical results 
obtained for a lattice with depth Vq = 9Ej., driven with frequency u/27i = 
1 kHz and scaled amplitude K/{hw) = 1.85. Under such conditions the single- 
band approximation is well justified; according to Eq. (!52|) . one has Jes = 
0.311 J. The drive was turned on quickly with a starting phase 0, as modeled 
by Eq. (jH]), and the images were taken 10 driving cycles later. Hence, provided 
that the solution (1M1) to the initial-value problem for the ideal gas discussed 
in Sec. O describes the real laboratory system correctly, Eqs. (17T]) and (1721) tell 
us that for any value of there should be stationary interference maxima at 
Piab = mod h{27i/d), or 

^ = mod 2 , (82) 

Pr 



using the recoil momentum pj- = hki^ as reference scale. This expectation is 
fully confirmed in Fig.O In all three cases considered there, with = 0, +vr/2, 
and — vr/2, the central peak is located at piab = 0, as corresponding to a con- 
densate occupying the Houston-Floquet state associated with the minimum of 
the quasienergy band (1751) . In addition there are side peaks at piab = ±2pi., cor- 
responding to the width h{2TT/d) = 2pr of the quasimomentum Brillouin zone. 
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Fig. 2. (a) Position of interference maxima observed in time-of- flight absorption 
imaging of condensates released from an optical lattice with depth Vq = 9 E,-, shaken 
with frequency uj/2tt = 1 kHz and scaled amplitude K/{huj) = 1.85, during one 
shaking cycle. The starting phases cp of the drive are (solid squares), +7r/2 (open 
squares), and —tt/2 (solid triangles); the images were taken 10 driving cycles after 
its turn-on. The maxima are centered around positions corresponding to momenta 
Piab = 0) ±2pr- (b) Modulation of the population difference of the side peaks during 
one shaking cycle. All three starting phases (symbols as in (a)) give rise to the same 
curve. 

The left panel of Fig. [2] depicts the time-resolved evolution of the positions of 
the respective interference maxima during a single cycle of the driving force. 
As anticipated, these positions remain practically constant in time, apart from 
apparent slight wigglings. The right panel of Fig. [2] shows the modulation of 
the height of the side peaks during one driving cycle, which stems from the 
oscillating time-dependence of the argument p of the factor |w(p)P determin- 
ing the envelope of the interference pattern in accordance with Eq. (!62|) . For 
all values of the starting phase we observe practically the same signals, which 
is in line with the surmise that they are signatures of the same state. 

Therefore, we interpret the data displayed in Fig. [2] as experimental signa- 
tures of Floquet condensates: The existence of a stationary, stable and lasting 
interference pattern indicates that the driven Bose gas actually tends to oc- 
cupy a single Houston-Floquet state. The fact that this pattern is centered 
around piab = 0mod2pr implies that this state is the one at the bottom 
of the quasienergy band ( 1751) . which, by virtue of Eq. ( 176|) . also is the one 
equipped with the lowest mean energy (fT3|) . This is not surprising, because we 
are starting with a condensate occupying the Bloch state k = 0, and because 
a drive of the form (HTj) does not change the quantum number k when it is 
turned on. 
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Fig. 3. (Color online) Bloch oscillations in driven optical lattices (Vb/i^r = 10), rec- 
ognizable as shifts of the interference pattern in time by one Brillouin zone of width 
2pj- within one Bloch cycle Tb. The Bloch frequency is set to wb/Stt = 242.4 Hz, 
whereas the driving frequency is ij/2iT = 3 kHz. The dimensionless driving ampli- 
tudes K/{hjj) correspond to 0, 1.5, and 4 (left to right). In (a), and to a lesser 
extent in (b), the pattern appears slightly smeared at times t/T^ = ±1/2, due to 
the occurrence of a dynamical instability. Note that the peak positions in (c) are 
centered around the edges itpr = zizhn/d of the quasimomentum Brillouin zone. 

If a Floquet condensate is to be regarded as an entity of its own, then it 
should behave as such when probed by a weak force acting on top of the time- 
periodic, driving one. This is actually the case, as illustrated by a sequence of 
further measurements summarized in Fig [3l Here we consider a lattice with 
depth Vq = 10 E^ driven by an oscillating force with an additional weak static 
component, written as 



m 



Fo- Ficos{ujt))e{t) 



(83) 



so that the time-dependent wave numbers fl69l) now become 

qk{t) = k ^ - r^ sm{ujt) 

h rko 



U) 



for t > 0, disregarding an initial phase from the outset. Therefore, the system 
undergoes a pure Bloch oscillation when Fi = 0: In response to a static force of 
strength Fq, one wave-number Brillouin zone of width 27^ /d is then traversed 
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at uniform "speed" Fq/U within the Bloch time 
2'nh 



Tr 



Fod ' 



^5) 



and the periodicity of the energy-quasimomentum relation (1701) in k gives 
rise to an oscilla ting •wave-pa cket motion, as long as interband transitions re- 
main negligible ( jZenerl . Il934[ l. We fix the Bloch frequency oj-b 1^2,71 = 1/Tb at 
242.4 Hz, and the driving frequency a;/27r at 3.0 kHz, more than 12 times 
higher. This clear separation of time scales implies transparent dynamics: 
Essentially, the oscillating component of the force f l83|) "dresses" the gas, 
making it condense into the Floquet state associated with the minimum of 
the quasienergy-quasimomentum relation fl7^ . as before. The relatively weak 
static component then merely probes this dressed system, making it behave 
as an undressed would if Eg. (1751) actually were the energy dispersion relation 



( lArlinghaus and Holthausl . 1201 ll ). 



The signatures of such dynamics are again visible in time-of-flight absorption 
images. We expect stationary interference peaks, corresponding to the min- 
ima of the dispersion relation ( I75p . while the argument of the envelope-giving 
function |w(j9)P in Eq. f l62p evolves in time as 



p_ 

Pr 



Plab 
Pr 



t 



K/{hw) 



vr 



sin((X't) 



^6) 



with K = Fid. This is precisely what is seen in Fig. |3l Here we show absorp- 
tion images obtained from condensates released at times ti/T-Q = i/12 with 
i = —6, —5, . . . , 5, 6, spanning one full Bloch cycle Tb (negative times here for- 
mally correspond to negative forces Fq in Eq. (!83l) ). In the leftmost panel we 
consider K/{tvjj) = 0, so that we a re dealing with the b are system undergoing 
an undisturbed Bloch oscillation ( jMorsch et al.l . 120011 ) . corresponding to an 
apparent displacement of the entire pattern by one Brillouin zone after one 
Bloch cycle. At t/T^ = ±1/2, when qkit) = ±7T/d and t he system passes th e 



20041 ) 



zone boundaries, a familiar dynamical instability occurs ( iZheng et al 
which becomes visible as a slight blurring of the peaks. 



The middle panel shows the corresponding pictures obtained for K/{hw) = 1.5, 
so that Jes = 0.512 J: Indeed this set of images looks quite similar to the pre- 
vious one, validating our view of a "dressed condensate" emerging through the 
application of the strong oscillating component of the force (I83p . and then re- 
sponding to the static component like a bare one with a renormalized hopping 
matrix element Jeff. The effective reduction of J by a factor of a bout one half 



also leads to a notable reduction of the dynamical instability (jZheng et al. 
2004h . 
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The third panel finally shows the images obtained for K/{huj) = 4, resulting 
in a negative effective hopping element, J^s = —0.397 J. This means that the 
minima of the quasienergy dispersion (1751) are now centered at the zone bound- 
aries k = ±71/(1, where the maxima of the original dispersion (1701) had been. 
Somewhat surprisingly, the gas spontaneously condenses into the associated 
Floquet state, leaving its signature in the form of two interference maxima 
neatly centered around p = ±pr- This strongly dressed system then again 
performs a Bloch oscillation, as wittnessed by the already familiar apparent 
displacement of the entire pattern by one quasimomentum Brillouin zone per 
Bloch cycle. 

Thus, our system actively seems to select the respective Floquet state with 
the lowest mean energy, and seems to conden se "by itself" into that state, as 



has also been reported by lLignier et al.l (120071 ) for the case of purely sinusoidal 



driving. This finding goes beyond the ideal-gas picture, and requires further 
studies. It is underlined by a related observation: If we instantaneously change 
the driving amplitude such that the absolute magnitude of Jes is preserved, 
but its sign is reversed, the original interference pattern caused by the initial 
state gradually vanishes, and the new pattern signaling the target state at the 
opposite band edge establishes itself within a few driving cycles. This process 
occurs in either direction, and appears to be fully reversible. This finding also 
implies that the dressed condensates settling down at the boundaries of the 
Brillouin zone when Jgs is negative actually are stable. 

Such experiments involving a sudden change of Jes are reminiscent of "quench- 
ing" experiments in which a parameter of a time-independent trapping poten- 
tial is varied instantaneously, and the ensuing rel axation dynamics are ob- 
served (JDziarmagal . l2010l : IPolkovnikov et al.l . 1201 ll ). However, it is uncertain 



at this point whether our observations fit into this picture, and the underlying 
relaxation mechanism — if it really is one — needs to be understood in detail, 
keeping in mind the fact that the system is isolated and thus cannot get rid of 
excess energy. Still, in view of the preliminary thoughts put forward in Sec. [21 
an interesting possibility suggests itself: The equilibrium state approached by 
a time-periodically driven, isolated Bose gas is determined not by usual ther - 



modynamics, but by "periodic thermodynamics" in the sense of lKohnI (120011 ) 
possibly involving a generalized temperature G which depends on — and can 
therefore, be controlled by — the external drive. 



7 Conclusions 



In a nutshell, the main results of the present study are encoded in Fig. [31 The 
observation that the middle panel of this figure closely resembles the leftmost 
one indicates that there are Floquet condensates, that is, macroscopically 
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occupied Floquet states of time-periodically forced Bose gases; the observation 
that the interference maxima are centered at the Brillouin zone edges when 
the effective hopping matrix element is negative, as seen in the rightmost 
panel, indicates that the Floquet condensate is carried by the Floquet state 
with the lowest mean energy. It remains to be seen whether these conclusions, 
drawn from one particular laboratory setting, also hold under more general 
circumstances. 

The observation of stable interference patterns in the presence of the time- 
periodic forcing gives strong support to the surmise, formulated in Sec. [21 that 
time-periodically driven, weakly interacting Bose gases confined by trapping 
potentals which forbid their escape effectively constitute equilibrium systems, 
rather than nonequilibrium ones; the equilibrium state being characterized by 
a constant distribution of Floquet-state occupation numbers. Under condi- 
tions such that the mean-energy constraint (IT^ is valid, the expected distri- 
bution is the Bose-Einstein distribution flTTj) . with the mean energies of the 
single-particle Floquet states replacing the single-particle energies appearing 
in time-independent situations. That distribution is characterized by two pa- 
rameters a and 7, corresponding to the generalized chemical potential z/ and 
to the generalized temperature 9 introduced in Eqs. fITS]) and flTIJl) . If this sug- 
gestion could be confirmed, it would open up further avenues of research: Also 
including Fermions, a time-periodically driven, isolated quantum gas should 
establish a temperature of its own, so that it might be interesting to explore, 
e.g., whether an already ultracold gas can be made even colder by applying an 
external time-periodic force. In any case, these tentative speculations clearly 
indicate that time-periodically driven quantum gases offer much more than 
mere visualizations of already known condensed- matter phenomena; in fact, 
they require the development of new concepts for time-dependent quantum 
many-body dynamics. 
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